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Evaluation of Exponentially Weighted Infinite Integrals by 
Quadrature (Hermite Quadrature) 
Closed subroutine, with one program parameter 
18 + 2N (see below) 

N(1«8 + T) milliseconds, where T is the duration in milli- 
seconds of the auxiliary subroutine. 
Location (may be used by auxiliary subroutine) 
When this routine is located at y, entry is made by the orders 



DESCRIPTION: 



p — aF 
^0 pF 

s 

p + 1 26 yF 
where a is the location of the auxiliary subroutine which 
computes the values of the function to be integrated., When 
control is returned to the right side of p + 1, the computed 
integral will be in the accumulator register and location 

y + 1^. 

To evaluate the integral 
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e"^ f (x) dx 



this routine uses a form of Gaussian Quadrature appropriate 

/ \ —x 

to the interval (- cd , od ) and the weighting function e : 



/ e""" f(x) dx -i^Z A^ f(x ) 
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(1) 



The values A, and x. are chosen^ in a manner such as to give udAJo 
k ■ . K " - 

truncation error when f (x) is a polynomialpOf degree 2N - 1 



or less. In the case where the factor e 
explicitly in the integrand. 
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does not occur 
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e"^ [e^ g(r)] dx ~ -^ g^^ A^^ e 
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2 K— ± p 

It is assumed that the function e g(x) may be closely 
approximated by a polynomial function » 

Because the actual values of the points x, and the weights 
A, and R may exceed 1^ they have been scaled down by powers 
of two- P and Q are defined in equations (l) and (2)_, and 
R is defined below in equation (3). 
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..[Note: Because the weights for these values of N become 
quite small, there may be appreciable roundoff errors for 
certain integrands ]« 

The auxiliary subroutine which computes f (x-, ) must take 
the scaling of these values of x, into account. The function 
values computed by the auxiliary are assumed to lie in the 
range -1 < f (x, ) < lo 



USE: 



SCALING: 



- 3 - 
The closed auxiliary subroutine is entered from the main 
routine with x, in the accumulator and link in Q| control 
is returned to the main routine with f (x, ) in the accumulator^ 
To use this routine the programmer copies the integration 
routine first on his program tape^ and immediately after 
it the parameters 5 points x , and weights A, or B, appropri- 
ate to his needs* These latter numbers appear on the tail 
of the library tape^ labeled by the number N of points at 
which the function is to be evaluated, and the type of weights 
(a, or B-, ) to be used« 

The scaling of the values of x, is such that the auxiliary 
subroutine is presented with x, * , where 



X, 



= 2"^^ Xj^, < R < 3, 



(3) 



y, ^ 

and the largest x, satisfied l/2 < (x, ) < 1 for all No 
k ' - ^ k 'max 

The computed integral is scaled down by 2 (or 2 }. 

For the convenience of the programmer^ the above scale factors 

are contained in the subroutine parameter at location y + l6, 

in the following form; 



ACCUEACY; 



(y+l6) = OR (y+l8) OP (y+l8+N), for A weights 

= OR (y+l8) DQ (y+l8+N) for B weights. 

The truncation error due to omission of powers of x higher 
than 2N is 



nl /jT" 4:,(2n) / \ , . = J- . / \ 
■ — -L__ f^ ■/ (zj^ where z is some point m (-qd , oojo 

2''(2n)l ; 

There will also be round-off errors which may be significant! 
for a discussion of these see the write-up of library routine 
E 5 - 195= 
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Preset point and weight addresses 



Plant subroutine entry and link 



clear svud. box 



Save LoSoP, of sma 



get X, I link to 



jump to auxiliary 



LoSoPo to A^ weight in Q 



X f (x, )_, accumulate M.SoPE 



step X, address 
step A, address 



exit via link 



temporary store for M.ByPo 
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Coefficient addresses and scale constants 
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end constant 



